Theory of heating of hot magnetized plasma by Alfven waves. 
Application for solar corona 
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The heating of magnetized plasma by propagation of Alfven waves is calculated as a function of 
the magnetic field spectral density. The results can be applied to evaluate the heating power of the 
solar corona at known data from satellites' magnetometers. This heating rate can be incorporated 
in global models for heating of the solar corona and creation of the solar wind. The final formula 
for the heating power is illustrated with a model spectral density of the magnetic field obtained by 
analysis of the Voyager 1 mission results. The influence of high frequency dissipative modes is also 
taken into account and it is concluded that for evaluation of the total coronal heating it is necessary 
to know the spectral density of the fluctuating component of the magnetic field up to the frequency 
of electron-proton collisions. 

PACS numbers: 52.35.Bj, 52.50.Sw, 96.60.P-, 94.30.cq 
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The discovery of the spectral lines of highly ionized 
iron in the spectrum of the solar corona showed in 1939 
that the corona is much hotter than the radiating surface 
of the Sun. Why it is hotter by about two orders of mag- 
nitude is still a very discussed problem, which survived 
not only the Second World War, but also the whole de- 
velopment of the modern astrophysics and satellite tech- 
nologies in the past half century. Definitely, this is the 
oldest still living mystery in front of the contemporary 
physics. 

On the other hand, heating of the solar corona is a phe- 
nomenon in the simplest physical system - a hot ideal gas 
of electrons and protons. All the phenomena in a com- 
pletely ionized hydrogen plasma are well-known and now 
it is a matter of a carefully proposed complementary set 
of experiments to arrive at the final conclusion. Accord- 
ing to the second law of thermodynamics each process 
in the plasma surrounding the Sun produces some en- 
tropy and can be considered as a heating mechanism of 
the solar corona. At least one new idea per year appears 
under the Sun, concerning the heating of its corona. In 
a review article on the subject one can find more than 
a hundred references. Often they are being classified as 
AC and DC models^ depending on whether a direct or 
an alternating electric current heats the plasma. But the 
coronal plasma can not be heated by an electric current, 
as it was proved thirty years age 8 : "...it is difficult to 
understand how a cooler gas (the electrons) can heat a 
hotter gas (the protons) within a region so far away from 
the Sun, which is the ultimate source of most of the en- 
ergy in interplanetary space" , see^. As far as the proton 
temperature is higher than the electron one T p > T e the 
heating mechanism acts on the protons and the viscose 
friction of the shear flow with a velocity field ~V(z,t) is 
a natural explanation for the bulk power Q of the dissi- 
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for divV = 0. Such a heating process can be observed 
in every turbulent fluid. Imagine a waterfall for which 
the amplitude of turbulent pulsations is smaller than the 
velocity of sound and as a result the related heating after 
the calming of the turbulent motion is smaller AT/T <C 
1 . That is why the turbulence in the solar plasma cannot 
give an augmentation by two orders of magnitude of the 
coronal temperature. For an effective heating mechanism 
we need an effective mechanism for the energy transport. 

This idea is supported by the recent development of 
the solar physics related to the random motion of mag- 
netic footprints. The creation of Alfven waves by the 
stochastic convection in the Sun is a very interesting hy- 
drodynamic problem for the complete heating scenario 
(for additional references see Ref. 3). The purpose of the 
present work is to analyze what kind of local satellites' 
measurements should be done in order to confirm or re- 
ject the Alfven wave mechanism for the coronal heating. 
Nowadays astrophysicists are able to use the available 
data provided by SOHO, TRACE, ACE, WIND in order 
to probe their model. However, the main defect of these 
experiments is the uncompleted set of data and that is 
why as far as the theory is concerned all the gladiators 
remain friendly sleeping in the arena. In the next sec- 
tion we will systematize the formulas for Alfven waves 
from well-known textbooks and in the final section we 
will analyze what kind of research has to be done. 



II. BASIC NOTIONS. MODEL EVALUATION 

Calculation of the dissipated power density is an in- 
dispensable ingredient in every theory which claims to 
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reveal the coronal heating mechanism. The real scien- 
tific problem is to explain the heating power, not the 
local temperature jump. The local temperature is only 
a parameter to be used in a realistic global model for 
the solar corona, which includes the acceleration of the 
solar wind by an expansion of the heated plasma and all 
processes for the local energy balance. 

We choose the z-axis oriented by the constant external 
magnetic field Bo = e z . We will investigate an Alfven 
wave, propagating on this axis. We can choose the y- 
axis to be perpendicular to the external magnetic field 
and the wave vector k = k x e x + k z e z . The Alfven waves 
are linearly polarized and the wave component of the 
magnetic field B = Bq + B is oriented along the y-axis 



B (r, t) = B e y , B = B yfi cos(k • r - u) A t)e~ 
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where the null index denotes the wave amplitude, loa 
is the frequency of the Alfven waves and the amplitude 
decrement 7 is supposed to be small enough. For the 
wave period averaged density of the magnetic energy we 
have 



We presume small fluctuations from the constant mag- 
netic field so that the maximal distortion angle of the 
magnetic field lines 
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is also small and we can use the linear theory of Alfven 
waves, propagating along the z-axis, whose phase velocity 
V A coincides with their group one v gr 



u A (k) = I V A ■ k| 
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V A k Zl k z > 0, vgr = V A - (9) 



For the bulk power density of the dissipated wave energy 
we obtain 
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In this model evaluation we have formally taken r = 0, as- 
suming that the weak space dependence of the magnetic 
field is included in the local value of the wave amplitude, 
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where /io = 4ir or 4ir * 10~ 7 in the practical system. 

If we follow the evolution of a wave-packet, the den- 
sity of the total energy (magnetic and kinetic) decreases 
exponentially over time 
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where the absorption rate of the Alfven waves 
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is given by a sum of the kinetic i^in and the magnetic 
i^mag viscosities, expressed in terms of the plasma viscos- 
ity ?y, density p and the electrical conductivity of the fluid 
(j; £o = 1/47T or Eq = I//10C 2 in the practical system. 

For an order of magnitude evaluation of the ampli- 
tude decrement^ we will consider only the Alfven waves 
loa = VAk z propagating along the external magnetic field 
lines k = fce z and discard the negligible influence of the 
resistivity is mag 
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where / = lo/2-k is the circular frequency. The relation 
between the velocity of the Alfven waves and the mag- 
netic field can be found in any textbook on magnetohy- 
drodynamics 



B' v (t)=B$cos(L0t), 
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measured by an AC satellite's magnetometer. Here the 
upper index (u>) emphasizes that up to now we have con- 
sidered a monochromatic wave whereas for the calcula- 
tion of the total absorption rate we have to perform a 
frequency summation over the whole Alfven waves spec- 
trum. In reality we ought to take into account also 
the influence of the Doppler effect, caused by the so- 
lar wind, but as a gedanken-experiment let us consider, 
that the frequency is measured in a system, bound to the 
coronal plasma. Imagine that the magnetometer has a 
frequency-filter and detects only frequencies in a narrow 
range (/i,/2), where f% — f% <C fi- The most impor- 
tant characteristic of the Alfven waves is their spectral 
density^, given by a sum of the wave energies in a narrow 
frequency interval 
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We can assume that this spectral density is given by an 
analogue filter, but in fact we have to do a Fast Fourier 
Transformation (FFT) for the magnetic field as a func- 
tion of time. Generally speaking <&(/) is a tensor, but 
for an illustrative purpose let us investigate only the y- 
components of the Alfven waves. Then the heating power 
is given by a frequency integral of the damping rate and 
the spectral density taken up to the maximal frequency 
for which the magnetohydrodynamic approach is still ap- 
plicable, i.e. the maximal electron-proton collisions fre- 
quency w max = l/r e : 
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In the damping rate we have to substitute the kinematic 
viscosity with that of a hydrogen plasma 6 
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is the Debye radius, T p and T e are respectively the proton 
and electron temperatures multiplied by the Boltzmann 
constant, and 



L„ = ln^»l 
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is the Coulomb logarithm as the thermal velocity of the 
protons is much smaller than the Bohr's velocity 



V T p < VBohr = -T, Mv Tp =T p . 
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Here we have used the elementary gas-kinetic evaluation 
for the viscosity 



?7 ~ n p MlpVT P , 
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as well as the mass density of the plasma, which is prac- 
tically determined by the density of the protons n p and 
their mass M 



p = Mn p . 



(19) 



In the derivation of the upper evaluation we have sup- 
posed that two protons with kinetic energy T p can be 
brought to a minimal distance r min 
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This distance parameterizes the effective cross-section <; p 
for the proton-proton collisions, which participates in the 
evaluation for the mean free path of the protons 



l p ^ p Tl p 1. 
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The Coulomb logarithm additionally increases the cross- 
section ? p 



?p ~ 7!"(2r m j n ) 2 Lp. 
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For the hot coronal plasma the kinematic viscosity is 
much bigger than the magnetic one i^kin 3> f ma g- This 
determines the diffusion of the magnetic field due to the 
ohmic resistance. 

For evaluating the maximal frequency at which mag- 
netohydrodynamics can still be applied we have to cal- 
culate the maximal frequency at which the density of the 
electric current, transferred by the electrons, j creates 



through the magnetic field B bulk Lorentz force, act- 
ing on the proton fluid with mass density p. This is in 
fact the frequency of the electron-proton collisions l/r e , 
which we take from the elementary kinetic gas theory. We 
can use the duration of the mean free path of electrons 
T e to determine the electric conductivity of the plasma 
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The mean time between two electron collisions is defined 
by the electron mean free path l e and their thermal ve- 
locity V T e 
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The relation between the electron density n e , their free 
path l e and effective cross-section <; e is the same as for 
protons 
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only the expression for the Coulomb logarithm should be 
changed 
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as even the velocity of the non-relativistic electrons is 
much greater than the Bohr velocity 
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The plasma frequency u p i is much higher than the colli- 
sions frequency 
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and yet it is located in the radio-range. 

These formulas give the well-known evaluations for the 
electron-proton collisions frequency 
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as well as for the conductivity of the fully ionized hydro- 
gen plasma 
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These model kinetic evaluations applied to the coronal 
plasma could be used for calculation of the heating power 
due to absorption of Alfven waves. For an analytical 
illustration we will utilize the interpolation formula, used 
for analysis^ of the magnetometric data from Voyager 1 
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Such a spectral density l// 2 is naturally explained in 
the framework of a stochastic Langevin MHD, as Mis- 
honov and Maneva pointed out 3 -, Eq. 41, and it is not 
related to any phase-coherent structures which could be 
distinguished from incoherent turbulent fluctuations, cf. 
with the work by Roberts and Goldstein^ 1 -. The Kol- 
mogorov turbulence f~ 5 ^ 3 give o 9 ' 10 a slightly (17%) dif- 
ferent power. The idea that convective granules produce 
a continuous stream of noise generating waves is actually 
very old 14 . We believe that this is a common property 
of all turbulent magnetized plasmas and we suggest the 
hypothesis that the 1/f 2 spectral density could be used 
even for the physics of accretion disks. 

Bearing in mind all above, the general formula Eq. (|13p 
gives the final model evaluation for the bulk heating 
power 
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This heating rate actually sets the upper limit for the 
absorbed bulk power that can be derived in the magne- 
tohydrodynamic framework when we take into account 
the influence of the diffusion Alfven modes with frequen- 
cies cu > l/r4 as well. However, these diffusion modes are 
able to heat only the narrow "skin layer" above the limb. 
Far away from the Sun (at the height of the satellites) we 
have to account only for those Alfven waves, propagating 
with lja > 1/ta. If we substitute the proton-electron col- 
lisions frequency cut-off 1 /r e in Eq. (fT3")) with the Alfven 
waves damping rate we will end up with a different as- 
sessment for the heating rate 
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Hereby we obtain another time constant related to kinet- 
ics of the protons that is rather significant for the theory 
of coronal plasma heating 
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This time constant Th together with the velocity of the 
solar wind v w - m d form a dimensionless parameter 



(35) 



which shows at what distance from the Sun rg un the dis- 
cussed mechanism for coronal heating stays important 
and should be incorporated in the coronal dynamics; at 
T <C 1 the influence of the Alfven waves heating is neg- 
ligible. According to our statistical approach, the heat- 
ing time-constant depends on the choice of the frequency 
cut-off. Therefore the role of the frequency cut-off for 
the coronal heating due to viscous damping of magneto- 
hydrodynamic waves is extremely important and we wish 
to find out what kind of physical processes in the plasma 
prevail to determine it. We have to pay attention that 



the ratio of these two cut-offs significantly varies with 
the solar radius, starting from a large number in the 
photosphere, where the collisions frequency dominates 
the Alfven waves attenuation, and reaching the order of 
unity in the lower corona, where the plasma (5 becomes 
rather small 
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As one can see, the most important characteristic in the 
result for the absorbed power Eq. ([32]) is the spectral 
density of the angle of deviation of the magnetic field 
lines from their state of static balance position. 
Due to this distortion as strings of a harp 
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we can listen to the impressive symphony of the convec- 
tive turbulence of the Sun and its corona is heated mainly 
by the high-frequency part of the spectral density. 



III. DISCUSSION 

A new theory often means data mining and experimen- 
tal data processing, but our purpose is not to extract the 
electron and proton temperatures from SOHO data; our 
purpose is to explain a local property of the magnetized 
plasma, pierced by Alfven waves. 

We have calculated the absorption power of the Alfven 
waves Qabs, that is the power transmitted per unit vol- 
ume from the Alfven waves to the protons. In order to 
incorporate this leading (according to our opinion) power 
in a coherent picture we ought to include Q a bs in the to- 
tal plasma energy balance which takes into account all 
possible energy losses or gains: the work due to the gas 
expansion, the change of the temperature which is pro- 
portional to the density of the kinetic energy, the en- 
ergy transmission from protons to electrons by collisions, 
bremsstralung of electrons, etc. In other words, Q abs 
should be the last missing ingredient of an almost com- 
pleted model for the coronal heating and creation of the 
solar wind. 

Few words should be said about how possibly our sce- 
nario could be rejected; the main defect of the contempo- 
rary coronal investigations is that none of the proposed 
mechanisms is consigned to the dustbin. Even the mag- 
netic reconnections which predict significant Ohmic heat- 
ing and T e » T p . As we propose Alfven waves as a 
heating mechanism the perpendicular to the static mag- 
netic field pulsations of the velocity and the magnetic 
field should be in relation corresponding to Alfven waves^ 



Eq. (69.11) 
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Such a relation was observed during the first space 
missions 12 . In other words the present heating scenario 
should be abandoned if the correlation coefficient be- 
tween the velocity and the magnetic field is not high 
enough. According to our picture the waves are gen- 
erated by the turbulence in the dense photospheric layer 
of the Sun where the influence of the magnetic field is rel- 
atively weak and mainly transversal incompressible mag- 
netohydrodynamic waves (MHD) waves are generated. 
These waves reach the regions with negligible pressure 
(i.e. low plasma (3 regions) but still their polarization 
corresponds to MHD waves in an incompressible fluid. 



IV. CONCLUSIONS 

The model problem for AlfVen waves transmission from 
a magnetized turbulent half-space to a non-turbulent 
half-space has been recently solved 3 . This model ap- 
proach can be further developed by an inclusion of the 
inhomogeneity of the magnetic field crossing the two half- 
spaces interface. However, for the coronal heating theory 
we consider as more appropriate to use experimental data 
for the Alfven waves spectral density. But the accuracy 
of the magnetic field measurements is much higher than 
that of the tachometry of the solar wind - that is why we 
focus on magnetometric data for the spectral density of 
the magnetic field Q(lo). This tensor ought to have two 
equal eigen values and another smaller one related to 
the static field component. Of course, it would be much 
better to have data for the magnetic field and velocity 
correlation. Investigations of the general 6x6 matrix, 
including all correlation coefficients, should be planned 
for the forthcoming space missions. The most impor- 
tant problem is where the spectral density of the MHD 
waves comes to an end. From a microscopic point of view 
the cut-off frequency is at electron-proton collisions rate 
1 /r e which is well in the range of the standard radio fre- 
quency measurements. However, the maximal frequency 
w max should be determined by the concrete physical con- 
ditions. For example, the inverse minimal size of the solar 
granulation could also determine the cut-off wave-vector 
k c . Another possibility for the maximal frequency is to be 
set by the waves propagation condition ojata > 1, which 
results in 

uj c = V A k c = vk 2 c = Va/v = B 2 /p i/p, k c = Va/v. 

(39) 

Above this frequency the waves are pure dissipative 
modes uj « iuk 2 with frequency-dependent depth of the 
"skin layer" S — \JvJuj. 

One of the problems we encounter today is that the 
high-frequency MHD waves are absorbed during the coro- 
nal heating and all the satellites' reports have a distorted 
spectral density like the red light of a sunset. In order 
to observe the bare MHD waves we need to have magne- 
tometric data from a satellite-kamikaze, but this would 
be a too expensive toy to be planned and realized. We 



need a more ingenious solution. We ought to analyze 
online which region of the Sun emits the MHD waves, 
investigated by the magnetometer. 

As a real break-through of the solar physics we consider 
the investigation of the correlation between the stochas- 
ticity of the space magnetic field delivered by the MHD 
waves and the stochasticity of the solar surface, obtained 
by optical data. Correlation of SOHO Dopplergrams 
with magnetograms is an important task, which should 
be put on the agenda in the global magnetohydrodynamic 
modeling of the solar corona. Before this is done we can 
just extrapolate the spectral density from the experimen- 
tally accessible region to the cut-off frequency l/r e , as 
it was done in the derivation of our model evaluation 
Eq. (fT3|) . This Q a bs can be compared with the necessary 
heating power of the global MHD models for the solar 
corona. This final formula contains only experimentally 
known parameters and can be used in the global MHD 
models for the solar corona. 

We would like to close with a qualitative discussion 
of the proposed mechanism in this work. As it is well- 
known, the leap of two orders of magnitude of the coronal 
temperature is rather sharp, compared to the height to 
the optically dense layers. This leap is a kind of a domain 
border, typical for the condensed matter physics. This 
sharp border ought to have at least qualitatively a good 
explanation in each mechanism for heating of the coronal 
plasma and this should be used as an important aesthet- 
ical criterion for the applicability of each coronal heating 
model. For most of them the closer we are to the dense 
layers the more intense is the heating and the observed 
temperature jump is practically unexplainable. Let us 
discuss this leap in the framework of the proposed sce- 
nario: slightly above the optically dense layers, emitting 
black body radiation tx T 4 , the dense plasma is transpar- 
ent for the MHD waves as the viscosity pv cx 77 cx T 5 / 2 is 
almost density independent. In this optically transparent 
layer the heating is negligible and the cooling is realized 
perhaps by a small residual convective turbulence. Imag- 
ine that all the plasma above is cool with T ~ 6 kK and 
let us consider its temperature stability. Imagine a local 
fluctuation - a small augmentation of the temperature. 
This increasing will lead to the enhancement of the vis- 
cosity r\ and the heat absorption from the MHD waves, 
piercing with an almost constant energy density through 
every layer of the plasma. This increasing of the heat- 
ing leads to an extra-augmentation of the temperature, 
which cannot be compensated form the thermal conduc- 
tivity in the rarefied plasma. Thus we have a positive 
feed back and the plasma explodes from the absorption 
of high-frequency MHD waves. This is simultaneously 
the mechanism explaining the creation of the solar wind. 
As a comparison let us mention that for none of the chem- 
ical explosions at room temperature we do not have an 
increasing of the temperature with as much as two orders 
of magnitude. It would be of a significant interest to in- 
vestigate the leap of the spectral density of the magnetic 
field at this frontier. 
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As a conclusion we would like to mention that in fact appreciated, 
the first historical explanation of the solar coronal heat- 
ing problem is still on the arena, retrieving another life 
or at least an impetus from the current satellites' data. 
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